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Flavor-Parity Breaking in the NJL Model with Wilson Fermions. 
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The Nambu-Jona-Lasinio model is chirally symmetric. Addition of a Wilson term explicitly breaks this sym- 
metry and leaves the model with a remaining parity-flavor symmetry. In the approximation of a large number 
of colors it has been shown that there is a phase where the remaining parity-flavor symmetry is spontaneously 
broken and that on the phase boundary all pions become massless. Using numerical simulations we conflrm the 
existence of this phase for the case of two flavors and two colors. We also conflrm the large N prediction of 
the existence of large flnite size effects that alter the shape of the phase boundary dramatically when periodic 
boundary conditions are used. 
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The lattice Nambu-Jona-Lasinio model is an 
interesting toy model of lattice QCD. For an ex- 
tensive analysis see [1] and references therein. 
The Lagrangian of the tvifo flavor model in 
Minkowski space and in continuum notation is: 
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The fermionic field '0/ is a fiavor SU{2) doublet 
and a color SU{N) vector, r = {ti,T2,T3} are 
the three isospin PauH matrices, ^ = 'y'^di^, and 
mo is the bare quark mass. The Lagrangian is di- 
agonal in color, in contrast with full QCD where 
it is diagonal in fiavor. The partition function is 
obtained after rotation to Euclidean space, addi- 
tion of a Wilson term and introduction of scalar 
a and pseudoscalar tt auxiliary fields (to make 
the action quadratic in the fermionic fields) [1]. 
The Wilson term breaks the chiral symmetry ex- 
plicitly down to fiavor symmetry. Therefore, if 
there are massless pions they can not be the Gold- 
stones of a spontaneously broken chiral symme- 
try. What is then the mechanism that provides 
massless pions? This question was posed in [2] 
for QCD and recently in [3] for the NJL model. 
There the NJL model is analyzed using large TV 
methods. The saddle point equations are: 
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with 



g{p, rrig) = ^ sin^ p^ + [rw{p) + rrig] 



w{p) =4-^cospi,, 
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where (Ts,Trs = \tTs\ are the uniform saddle point 
fields, /3i = j^, r is the coefficient of the Wilson 
term (we will set r = 1), and rriq is the constituent 
quark mass. There are two solutions to these 
equations corresponding to two phases. (I) A 
fiavor-parity symmetric phase where the order pa- 
rameter TTs = {'^ ('^«75T^)s = 0) and massive 
pions m„j^,m„2,m„^ ^ 0. (II) A spontaneously 
broken fiavor-parity symmetry phase where, for 
example, tt^ = (0,0, tt^), tt^ ^ (-^ (*«7bT*)s ^ 
0), two pions are massless m„^ = m„^ = 0, and 
one is massive m^g ^ 0. The two massless pions 
are the Goldstones of the spontaneously broken 
fiavor symmetry. The mass of the third pion can 
be interpreted as the inverse correlation length 
associated with the spontaneous breaking of par- 
ity symmetry. On the phase boundary this dis- 
crete symmetry undergoes a second order Ising- 
Hke phase transition where the associated corre- 
lation length diverges and therefore the mass of 
the third pion vanishes. So the phase boundary, 



which is a line in the /3i , toq plane, is characterized 



by TTs 



0, niT, 



0. This phase 



diagram, determined from the saddle point equa- 
tions, was calculated in [3] for an infinite volume 
lattice using a Newton procedure to perform the 
integrals. Here we will concentrate our efforts in 
finite size lattices for two reasons. First, we will 
perform numerical simulations for finite TV = 2 
in order to confirm the existence of the sponta- 
neously broken phase. Second, on finite lattices 
the zero fermionic modes alter the shape of the 
phase diagram significantly. We study this using 
large N methods and numerical simulations. 

On a finite lattice the phase diagram can be 
computed from the saddle point equations at the 
transition point (tt^ = 0): 



\" 



1 



mo 



LlLt ^ gik,mg] 



+ 



\^ 



w{k) 



Pi 



2N 



mo 



LIL 



LlLt ^ gik,mg) 
y^ mg + rw{p) ^p 



t— 9iP,mg) 



(4) 



(5) 



For given mq we calculate mo from 4. Insert- 
ing this into 5 we find /3i. On a periodic lat- 
tice the momentum takes values pi = 2j2i ^ fi- = 
0, 1, • • • Lj — 1. Then g{p, nig) vanishes at: 
mg= 0, p= (0,0,0,0), 
iriq = —2, p = (tt, 0,0,0) and permutations, 
mq = —4, p = (7r,7r, 0,0) and permutations, 
rUg = —6, p = (7r,7r,7r,0) and permutations, 
TUg = -8, p = (7r,7r,7r,7r) . 

The 16 momenta where g{p,mg) vanishes corre- 
spond to the origins of the Brillouin zones of the 
16 doubler species. As a result, at these values of 
mq, the corresponding terms of the momentum 
sum in eq. 5 are singular, and therefore so is the 
function /3i(mo) that describes the phase bound- 
ary. This is shown for an 8'^ x 16 size lattice by 
the solid line in Figure 1. 

Inside the octopus-like graph the fiavor-parity 
symmetry is spontaneously broken. At point "1" 
of this figure m,, is very large and positive. As 
m,g decreases, crosses zero and tends to very large 
negative values, we transverse the whole solid line 
from "1" to "2" to "3" • • • to "12" where m, is 
very large and negative. The "prongs" of this fig- 
ure, points "3", "5", "7", "9", and "11", extend 



all the way to /3i = oo where they correspond to 
m,q = 0, —2, —4, —6, and — 8, respectively. When 
the singular terms described above are neglected 
from the calculation (the contribution of these 
terms disappears in the infinite volume limit) we 
obtain the dotted line of Figure 1. This line is 
quantitatively very close to the infinite volume 
phase line of [3] . As the volume increases to infin- 
ity the "spikes" of the solid line shrink in "thick- 
ness" and eventually are mapped close to the dot- 
ted line, producing the infinite volume phase line 
of [3]. These spikes are a very strong finite size 
effect. For example, at finite volume, in order to 
obtain a very small niq, one has to be near point 
"3", i.e. at very large /3i. As the volume becomes 
infinite very small values of nig can be obtained 
at finite /3i close to point "2". For more details 
and figures see the second reference in [1]. 




Figure 1. Phase diagram with periodic BC 

When antiperiodic BC are used even along one 
direction (say along the time direction) g{p,mq) 
has no zeros (for example the smallest momentum 
is (0,0, 0, -£-)). As a result, /3i(mo) has no singu- 
larities. The spikes are still present but they do 
not extend to /3i = oo. Their extension along the 
/3i direction is roughly proportional to j^ ^\ ■ 
The phase diagram for lattice sizes 16^ x 32, 
16^ X 64, 8^ X 16, 8^ x 64 is shown in Figure 2. In 
the order mentioned above the 16^ x 32 has the 
smallest size "spikes" and the 8^ x 64 the largest. 
In order to confirm the large N predictions we 
have performed Hybrid Monte Carlo simulations 
on the CM-2 supercomputer at SCRI on 8"^ x 16 



lattices at A'' = 2. In the parity broken phase we 
define the "magnetization" 



M„ = '^TT^^t, T^x,t = i^x,t'y5T^x,. 
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The massive pion propagator is measured from 
the decay in the "t" direction of the correlation 
function c(t) of the component of tt along M„. 
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Figure 2. Phase diagram with antiperiodic BC 
along the time direction 

The results for the pion mass at /3i = 2.5 are 
denoted by crosses in Figure 3. The solid line is 
the large N prediction. This confirms both the 
existence of a spontaneously broken flavor-parity 
phase as well as the presence of a "spike" at that 
region of the phase diagram. The same plot, but 
for /3i = 0.5, is shown in Figure 4. In this case, 
as can be seen from Figure 1, when mo is varied 
from —3 to —1.5 one crosses the phase line three 
times. The crossing at mo « —2.76 is predicted 
fairly well by the numerical results. The other 
two crossings entail a very small pion mass and 
are missed by the numerical results (albeit some 
indication from the shape of the curve) . Similar 
conclusions can be drawn by looking at the plot 
of |M„| = M„ in Figure 5. 

In conclusion, both large N predictions, namely 
the existence of a spontaneously broken flavor- 
parity phase as well as the presence of the 
"spikes" on flnite lattices with periodic boundary 
conditions, were verifled numerically at TV = 2. 
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Figure 3. /3i = 2.5, lattice size is 8'^ x 16 
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Figure 4. /3i = 0.5, lattice size is 8'^ x 16 
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Figure 5. j3i = 0.5, lattice size is 8^ x 16 
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